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The origin of primary nitrogen in galaxies 
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Abstract. We investigate the role of stellar axial rotation on the nitrogen nucleosynthesis at low metallicities Z. 
For this purpose, we have calculated models with initial masses between 2 and 60 Mq at Z=0. 00001 from the 
zero age sequence to the phase of thermal pulses for models below or equal to 7 M0 , and up to the end of central 
C-burning for the more massive stars. The models include all the main physical effects of rotation. We show that 
intermediate mass stars with rotation naturally reproduce the occurrence and amount of primary nitrogen in the 
early star generations in the Universe. We identify two reasons why rotating models at low Z produce primary 
^'*N: 1) Since the stars lose less angular momentum, they rotate faster. Simultaneously, they are more compact, 
thus differential rotation and shear mixing are stronger. 2) The H-burning shell has a much higher temperature 
and is thus closer to the core, which favours mixing between the two. 
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1. The context: observations and models 

For more than 20 years, the nucleosynthetic origin of pri- 
mary nitrogen, which is one of the most abundant and im- 
portant elements in the Universe, for life in particular, has 
remained a deep mystery in astrophysics (Edmunds and 
Pagel |1978| Barbuy |l983t Ca rbon et al. |l987t T huan et al. 
|1995| ; Izotov & Thuan [1999| ; Henry et al. ^OOOD . Nitrogen 
is mainly produced by the CN cycle of the CNO reactions 
which catalyze hydrogen burning in stars. In the CN cy- 
cle, the reaction -'^'*N(p,7)^^0 which destroys nitrogen has 
a very small cross section, which allows "'^^N to accumu- 
late with time (e.g. Arnett 1996). Thus, ^''N is usually the 



daughter element, hence a secondary element, of the CNO 
initially present in stars. If nitrogen is secondary, the in- 
crease in the abundance of ^''N should be proportional to 
the initial CNO content and consequently the ^^N-content 
will be proportional to the square of the CNO and metal 
content in a galaxy. If ^"^N comes from a primary pro- 
duction, i.e. is produced from the original hydrogen and 
helium, the ^''N-abundance is proportional to th at of the 
other primary heavy elements (Talbot & Arnett 1974 ). 

The first evidence of a primary production of ^"'N in 
the early phases of the evolution of our Galaxy came from 
the study of the very old an d low metallic ity sta rs. Several 
authors ( Edmu nds & Pag el |l978| ; Barbuy [l983l Tom kin & 
Lamb ert |1984| ; Matteucci |l986| ; Carbon et al. |l987t Henry 
et al. ^000 ) have shown that the ratio N/O of nitrogen 
to oxygen remains constant with a plateau at log N/O 



~ -1.6 in the early evolution of the Galaxy, thus imply- 
ing a primary origin of nitrogen. Later when the metal 
abundance has reached about 1/5 of the solar abundance, 
the N/C and N/O ratios grow rapidly, as expected for 
secondary elements. Thus, primary production seems to 
occur only in the very early phases of galactic evolution, 
while secondary production appears later, when the metal 
content of a galaxy is higher. A second piece of compelling 
evidence is provided by the study of the N/O ratios in ion- 
ized HII regions of blue compact dwarf galaxies (Thuan 
et al. |1995| ; Izotov & Thuan |1999| ; Izotov & Thuan |00^. 
There also N/O is observed to be constant at very low 
metallicities Z. The blue compact galaxies are in general 
galaxies where the average star formation rate has been 
very low in the past and which are therefore still in an 
early stage of their chemical evolution. An example is the 
galaxy IZw 18, which has the lowest known metallicity 
(1/50 of solar), and which shows indications of primary 
nitrogen (Kunth et al. |1995| ; Izotov & Thuan [19991 ). 
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A third indication for the presence of primary nitrogen 
comes from the observed N/O gradient in spiral galaxies. 
If nitrogen is a purely secondary element, the N/O gradi- 
ent should be equal to the gradient of heavy elements, and 
in particular to the gradient of O /H in these galaxies. The 
fact that the N/O gradient is relatively flat at low metal- 
licity Z indicates that production of nitrogen is dominated 
by primary processes a t low Z and secondary proc esses at 
high Z (Ga rnett et al. |l997t Ferguson et al. |l99j ; Henry 
& Worthey 1999 ). A problem has arisen some years ago, 
because some low Z damped Lya systems, which are likely 
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protogalactic structures at relatively high redshifts, have 
N/0 ratios lower than those observed in the HII regions 
of blue compact dwarf galaxies of the same Z (Pettini et 
al. 1995). However, the apparent discrepancy has been re- 
solved by models of damped Lya systems which account 
for both ionized and neutral regions (Izotov et al. 2001). 

Observationally, it is still uncertain whether the ob- 
served primary nitrogen originates from massive or from 
intermediate mass stars. The claims that primary nitro- 
gen is made in massive stars are mainly based on the low 
scatt er of the obse rved N/0 ratios at lo w Z ( Matteucci 
1986| ; Thuan et al. |l995t Izotov & Thuan |l999| ; Izotov & 



Thuan pOOOl) . The argument is that, if nitrogen is synthe- 
sized in massive stars, there is no time delay between the 
injection of nitrogen and oxygen and thus a rather small 
scatter would result. On the contrary, if the primary ni- 
trogen is made in intermediate mass stars, the N/0 ratio 
increases with time, since these stars release their nitro- 
gen much later than the oxygen made in massive stars. 
This would lead to a larger scatter in the observations of 
galaxies at various stages of their evolution. Some recent 
studies found that a significant scatter does exist (Garnett 
19901 ; Skillman et al. |1997D . Also, Henry et al. ( |2000| ) have 



calculated models which support the view that intermedi- 
ate mass stars between 4 and 8 M0, with an age of about 
250 Myr, are likely to dominate the nitrogen production. 

Let us examine what stellar models tell us. The con- 
ditions needed for the production of primary nitrogen are 
very simple. In a star which has both a helium burning 
core and a hydrogen burning shell, some amount of the 
new carbon synthesized in the core must be transported 
into the hydrogen burning shell, where the CNO cycle 
will convert it into primary ^^N . For massive stars, only 
models with ad hoc hypotheses are able t o do this. For 
example, in some models (Timmes et al. 1995 ), mixing 
has been arbitrarily "permitted" between the helium- and 
hydrogen-burning zones. Without any physical explana- 
tion, it is difficult to understand why the production of 
primary nitrogen only occurs at low metalliciti es. M odels 
of metal free Population III stars (Umeda et al. 200C| ) pro- 
duce some primary nitrogen, but in too low quantities. Up 
to the phase of thermal pulses, asymptotic giant branch 
(AGB) stars produce no primary nitrogen. Only during 
thermal pulses, some helium products might be trans- 
ported into the hydrogen burning shell to produce primary 
nitrogen. However, the existing AGB models (Renzini & 
Voli 1981; Marigo 1998) are "synthetic" models, which 
means that the model parameters follow some analytical 
relations that have been fitted to the observations. While 
this may be useful for some purposes, it cannot be claimed 
that it represents consistent physics leading to primary ni- 
trogen production. 

2. Primary nitrogen production in rotating models 

We construct stellar models with axial rotation and very 
low metallicity Z, typical of the early star generations. 
Rotation is clearly important, because there are indica- 




Fig. 1. Evolution of the distribution of angular velocity Vt 
as a function of the distance to the center in a 20 M© star 
with an average velocity of 230 km s~^ and Z=0. 00001. 
Xc is the hydrogen mass fraction at the center. The dotted 
line shows the profile when Xc=0. The dashed line shows 
the same at Z=0.004 when ^^=0.28. 



tions from studies in the Magellanic Clouds that stars 
ma y on the average rotate faster at lower Z (Maeder et 
al. 1999). Our models include the main effects of rotation 
on stellar evolution: the hydrostatic effects, the merid- 
ional circulation which transports both angular momen- 
tum and chemical elements (Maeder & Zahn |1998[ ), the 
effects of horizontal turbulence (Zahn 1992| ), the effects of 
shear mixing (Zahn [1992| ; Maeder & Meynet |2001| ) and 
also the effects of rotation on the mass loss rates (Maeder 
& Meynet 2000] ). Such model physics has already been ap- 
plied to stars of solar composition and to the Magellanic 
Clouds. In both cases these models have resolved major 
discrepancies (Meynet & Maeder 2000; Maeder & Meynet 

pooi| ). 

We have calculated models of 2, 3, 5, 7, 9, 15, 20, 40 
and 60 M© at Z =0.00001 from the zero age sequence to 
the phase of thermal pulses for models below or equal to 7 
Mq , and up to the end of central C-burning for the more 
massive stars. Fig. 1 shows an example of the evolution 
of angular velocity Vl until the end of the Main Sequence 
phase. The interesting property of very low Z models is 
that the 17-gradients are much steeper than at higher Z, 
because of the lower mass loss rates at lower Z and the 
greater compactness of the star. The consequence of the 
larger f2-gradients is that there is more mixing due to 
the shear instability and thus shallower internal molec- 
ular weight gradients (/i-gradients). Comparisons of the 
surface enrichment in A-supergiants in the SMC and in 
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Fig. 2. Distribution of the main elements in the region 
of the hehum- and hydrogen-burning sheUs during the 
TP-AGB phase of an initial non-rotating 3 Mq star at 
Z^O.OOOOl. 
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Fig. 3. Same as Fig. g for a 3 Mq rotating model dur- 
ing the TP~AGB phase. The average rotational velocity 
during the Main Sequence is 230 km s~^. 



burning shell. Later, on the AGB, the zone between the 
He- and the H-burning shells, enriched in primary nitro- 
gen, shrinks, but at the same time more and more primary 
is entering the outer convective zone. This primary 
^^N in the outer envelope will escape from the star, either 
in the stellar wind or during the formation of a plane- 
tary nebula. This mode of nucleosynthesis is quite different 
from the classical scenarios in which the primary nitrogen 
is produced during the thermal pulses (TP) AGB phase. 
Fig. 2 and 3 show the abundances of some nuclei in the 
intershell region in TP-AGB models of 3 Mq models, with 
and without rotation, at the luminosity of the 5th thermal 
pulse. We notice first that the CO core is slightly larger 
with rotation and that the chemical gradients at its outer 
edge are also shallower. The most remarkable difference is 
the much higher "'^'^N-content in the intershell region (i.e. 
where helium has its highest plateau) and in the adjacent 
external deep convective zone. This ^"^N-content is orders 
of magnitude higher than the total abundance of original 
CNO elements and it is thus of primary origin. The models 
between 2 and 7 Mq show a similar behaviour. We em- 
phasize also that similar models with rotation at Z=0.020 
do not produce any primary nitrogen (Meynet & Maeder 
[2000 ). The mechanism which produces primary nitrogen 
only works efficiently for Z less than ~ 1/5 of solar. 

There are two physical reasons for this production of 
primary "'^^N in rotating models at low Z . 1) As explained 
above, the r2-gradients are steeper, while the /x-gradients 
are shallower at lower Z. These effects are present during 
the whole evolution and enhance shear mixing. Thus, at 
very low Z with rotation, there is a very efficient transport 
of new ^^C from the core to the H-burning shell, where the 
CNO cycle turns it to primary ^^N. 2) The second reason 
is that at lower Z the CNO burning occurs at much higher 
temperature T. Indeed, the amount of energy produced in 
TP-AGB stars is of the same order of magnitude, what- 
ever the value of Z. In order to maintain the rate of energy 
production when, as in the present models, the CNO con- 
tent is about 2 • 10^ times smaller than solar, T has to be 
about (2 • 10^)'-'^/^) — 2.3 times higher, because the energy 
production of CNO is proportional to at the typical 
T = 8 • 10"^ K, where H-shell burning occurs in TP-AGB 
stars of very low Z. Since the H-burning shell needs to 
be at much higher T in lower Z stars, it is much closer 
to the edge of the He-burning shell, which is sitting on 
the outer edge of the CO core and which produces ^^C. 
The value of T in the CO-core is not very sensitive to the 
initial composition. Thus, the intershell distance is much 
shorter at very low Z (by about an order of magnitude in 
the present case) and the transport of ^^C from one shell 
to the other is easier. 



the Milky Way fully support the presence of more mixing 
at lower Z (Venn |1999| ; Maeder & Meynet pOOl] ). 

The combination of the effects of low Z and rotation 
enables primary nitrogen production during the core He- 
burning phase, in the zone between the core and the H- 



3. The nitrogen stellar yields at low metallicity 

Fig. 4 shows the amount of nitrogen, the so-called stellar 
yields mpNi4 as a function of the initial stellar masses. We 
notice that intermediate stellar masses of 5 to 7 Mq have 
the largest stellar yields of primary nitrogen. When these 
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Fig. 4. Amount of nitrogen (mainly primary in rotating 
models) produced as a function of the initial stellar mass. 
The quantities given here are the so-called stellar yields 
mpNi4, where m is the initial mass and pni4 its fraction 
in form of total ejected "'^^N. The initial rotation velocity 
Vini on the ZAMS is indicated. For the initial mass stars 
between 3 and 15 M©, the average rotational velocity dur- 
ing the Main Sequence is around 230 km s~^. The yields 
weighted by the Salpeter mass spectrum is shown by a 
dotted line (arbitrary units). 



yields are weighted by the Salpeter IMF (Fig. 4), we see 
that the main contributors are between 2 and 5 Mq. For 
the initial stages of the evolution of galaxies considered 
here, we may to first approximation take the ratio of the 
net yields in nitrogen and oxygen as an estimate for the 
observed N/0 ratio. For the production of oxygen, we use 
the relation between the mass Ma of the He-C-0 core 
at the end of the He or C- burning phase and the yield 
in oxygen by Arnett ( 1991 ). Taking both the N and O 
productions for the models at Z = 0.00001 and convoluting 
them with a standard initial mass function, we get a ratio 
log N/0 = -1.88. This result is within a factor of two of 
the value log N/0 ~ -1.6 observed at low Z. This is a 
very satisfactory agreement for two reasons. 1). Here, we 
have considered models with an average velocity of 230 
km s~^. However, the production of primary ^"^N increases 
rapidly with rotation, thus the production for the average 
velocity is smaller than the average production for the 
actual distribution of rotational velocities. Moreover, there 
seem s to be more fast rotators at lower Z (Maeder et al. 



1999|). 2). Nitrogen is ejected mainly by AGB stars with 
ejection velocities of a few 100 km s^^, while oxygen is 
ejected by supernovae at much higher velocities of 10'* km 
s~^ or more. Thus, a fraction of the oxygen produced may 



escape from the parent galaxy, leading to a higher N/0 
ratio than in the simple estimate made here. 

4. Conclusion 

For the first time, we present models of star evolution with 
rotation and very low metallicities, typical of the first stel- 
lar generations in the Universe. We show that stellar ro- 
tation in very low Z models leads to the production of 
primary nitrogen in amounts that are in global agreement 
with those observed. The main sources are intermediate 
mass stars with initial masses between 2 and 5 M© (see 
Fig. 4). The yields of s-process elements at low metal- 
licity will also be deeply affected by rotation (Langer et 
al. 199£). Indeed, as can be seen from Figs. 2 and 3, the 
abundance of the main neutron source ^^C is considerably 
enhanced in the intershell region of the rotating model. 
Finally, we emphasize that a certain amount of primary 
^^Ne should also be synthesized as a result of the burning 
of primary ^*N in the He-burning shell. However, con- 
trary to primary ^''N, the exact yields in ^^Ne are difficult 
to predict, since they critically depend on the mass loss 
and mass cut leading to the final remnants. 
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